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1 Introduction 



The article gives an approach to the exponential stability of equations of the form 

8 V + AU = F, (1) 

where we denote the derivative with respect to the temporal variable by do. The (unbounded) 
linear operator A, acting on some Hilbert space, is assumed to be maximal monotone and 
we refer to the monographs [?| [15] for the topic of monotone operators on Hilbert spaces. 
Equation ([1]) is completed by a constitutive relation of the form 

V = MU, 

where M is a bounded linear operator acting in time and space. Thus, we end up with an 
equation of the form 

(d M + A)U = F, 

which we refer to as an evolutionary equation. This class of problems was introduced by 
Picard in |16J , where A was assumed to be skew-selfadjoint, and it was illustrated that many 
equations of classical mathematical physics are covered by this abstract class. The well- 
posedness for such problems is proved, by showing that the operator <9o-M + A is boundedly 
invertible in a suitable Hilbert space. More precisely, the derivative do is established as a 
normal, continuously invertible operator in an exponentially weighted L2-space and M is 
defined as a function of d^ 1 (see Section 2) and the well-posedness is shown under a positive 
definiteness constraint on the operator M (see Solution Theory] and Theorem 12.41 of this 
article). Moreover, the question of causality, which can be seen as a characterizing property for 
evolutionary processes, was addressed, which leads to additional constraints on the operator 
Ai, namely that M is defined via the Laplace-transform of an analytic and bounded function 
M : Bc(r,r) — > L(H) (for more details see Section 2). Especially the analyticity of M is 
crucial for the causality, due to the correlation of supports of /^-functions and the analyticity 
of their Laplace-transforms by the Paley- Wiener Theorem (see [20^ Theorem 19.2]). Later on 
these results were generalized to the case of A being a maximal monotone relation in [21] [23] . 

In this work we give sufficient criteria in terms of the function M, which yield the exponential 
stability of the evolutionary problem. The study of stability issues for differential equations, 
which goes back to Lyapunov (see |14| for a survey) , has become a very active field of research 
for many decades and there exist numerous works dealing with this topic. We just like to 
mention some standard approaches to exponential stability. The first strategy goes back to 
Lyapunov. The aim is to find a suitable function (a so-called Lyapunov function) yielding 
a certain differential inequality, which allows to derive statements about the asymptotic be- 
havior of solutions of the differential equation. A second approach, which applies to linear 
differential equations, is based on the theory of semigroups. In this framework different cri- 
teria for exponential stability were derived in terms of the semigroup or its generator, e.g. 
Datko's Lemma ([8] or [9] p. 300]), Gearharts Theorem ( [11] or [9] p. 302]) or the Spectral 
Mapping Theorem (see p. 302, Theorem 1.10]). A third approach uses the Fourier- or 
the Laplace-transform of a solution to derive statements of their asymptotics. These meth- 
ods are sometimes referred to as Frequency Domain Methods. In our framework it seems 
to be appropriate to employ the last approach, since, by the definition of M, methods of 
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2 The framework of evolutionary equations 

vector-valued complex analysis are at hand through the Laplace-transformation. Note that, 
due to the general structure of evolutionary equations, the results apply to a broad class of 
differential equations, such as differential-algebraic equations, equations with memory effects 
or integro-differential equations, where semigroup methods may be difficult to apply. 

The article is structured as follows. In Section 2 we recall the framework of evolutionary 
equations and its solution theory (Theorem I2.4p . Section 3 provides an abstract condition for 
exponential stability for evolutionary equations in terms of the function M and the proof of 
the main stability result (Theorem I3.T() . To illustrate the versatility of the previous results, 
we discuss several examples in Section 4. We begin to study differential-algebraic equations 
of mixed type and derive a condition for their exponential stability. Moreover, we give a 
concrete example for such an equation, which seems hard to be tackled by other approaches, 
since the type of the differential equation switches on different parts of the underlying do- 
main. Furthermore, we give a possible approach how to deal with initial value problems. In 
Subsection 4.2 we consider an example of a partial differential equation with memory effect. 
A similar problem was also treated by Batkai and Piazerra in [2], using a semigroup approach. 
In contrast to their result our approach directly extends to the case of differential-algebraic 
equations with delay. We conclude the article by discussing a parabolic integro-differential 
equation with an operator-valued kernel. We adopt the ideas of |22j to derive the exponential 
stability. 

Throughout let H be a complex Hilbert space. We denote its inner product by (-|-), which is 
assumed to be linear in the second and conjugate linear in the first argument. We denote its 
induced norm by | • |. 



2 The framework of evolutionary equations 

We recall some basic notions and results on linear evolutionary equations, i.e. equations of 
the form 

(d , e M(d^) + A)u = f, (2) 

where A : D(A) C H — > H is a linear, maximal monotone operator, <9o,p denotes the time- 
derivative, established in a suitable Hilbert space and M(8q ) is a bounded linear operator 
in time and space, a so-called linear material law. We refer to \16\ W7\ [T3"l [2T| [23] for more 
details and proofs of the following statements. First we begin to introduce the Hilbert space 
setting, where we want to consider equation ([2]). 

For pGlwe define the space H Qt o(M.; H) as the space of (equivalence classes of) measurable 
function / : R — > H which are square-integrable with respect to the exponentially weighted 
Lebesgue-measure e~ 20t dt, equipped with the inner product 

(f\9)H eA R;H) ■■= J (m\g(t))e-^ dt. 

E 

Note that Ho t o(R; H) is just the space L2(M.;H). We define the derivative <9o >£1 as the closure 
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of the operator 

do,g\cg°(R;H) '■ C^°(R;H) C H Bi q(R;H) — > H gi0 (R;H) 

4> 0', 



where we denote by C£°(M;H) the space of function : R — > H, which are arbitrarily often 
differentiate and have compact support. 

Remark 2.1. 



(a) For each g£l the operator do, s is normal with *Rt do >e = ^ [Qq^ + <9q g j = g. Moreover we 
obtain that <9q = — <9o, e + 2g. In particular, for £> = the operator <9o,o is skew-selfadjoint 
and coincides with the usual weak derivative on L2(JS.;H) with domain Hoi(M.;H) = 
H^R-H) = W%(R;H). 

(b) For g ^ the operators <9o,g on H g o(R; i?) and <9o,o + £ on ^(Rj are unitarily equivalent 
and c?o,£ is boundedly invertible with = t^t. For £> = the operator c?o,o is n °t 
boundedly invertible, since is in its continuous spectrum. 

(c) For g 7^ the inverse operator is given by 

/loo u ( s ) ds if £ > °> 
— u(s) ds if q < 



(5>) (o 



for u G i/ e) o(R; ^0 an d almost every t G R. This representation yields, that for g > the 
operator c^T* is forward causal whereas it is backward causal for g < 0. 

(d) Let £ G R and denote by C Q : i7 ei o(M;i?) — > L2(M>;H) the so-called Fourier- Laplace 
transformation, defined as the unitary extension of the operator given by 



(C e <p) (*) = -L J e-^e-^t) dt 



for <p G C C °°(M; #) and x G R. Then 



where by m : D{m) C L2(R;H) — > L2(R;H) we denote the multiplication-by-the- 
argument operator with maximal domain, i.e., (mf) (t) = tf(t) for almost every t € M. 
and every / G D(m) := {/ G L 2 (R; il) | (i H> i/(i)) G L 2 (R; i?) } . Note that in the case 
Q = 0, ([3]) is just the usual spectral representation for the weak derivative on L2(M.;H) 
via the Fourier-transformation (see [H p. 112]). 

We consider ([2]) as an equation in the Hilbert space H gi o(M.; H). As a matter of physical 
relevance, we require that the corresponding solution operator fdo jS M(dQ^) + A^j , if it 

1 A mapping F : D(F) C i? e;0 (R; H) — > J/ e ,o(R; ff) is called forward causal, if for each ■D(.F) with / = g 

on some interval ] — oo, a[ for a £ R the functions F(f) and -FXg) coincide on the same interval ] — oo, a[. 
Analogously F is called backward causal, if for each f,g G D(F) with f = g on some interval ]a, oo[ for 
o£R the functions and F(g) coincide on the same interval ]a, oo[. 
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2 The framework of evolutionary equations 



exists, is forward causal, that means, roughly speaking, that the present state of the solution 
u does not depend on the future behavior of the source term /. With Remark 12. II (c) in mind, 
we therefore assume that g > 0. 

We now define the operator M(<9^) for g > with the help of formula ([3]). 
Definition 2.2. Let r > and M : B c (r,r) ->■ L(H) H For g > ± we define 

M ( a-):=0/(-L_)/: e , 

where M is defined as the multiplication operator (m ( r^f^ ) f\ (t) '■= M (j^^j fit) 

on L 2 (R; H) with domain {/ G L 2 (R; H) (t M /(*)) G L 2 (M; H) }. We call M(d^) 

a linear material law if the function M belongs to 7i°°(Bc(r,r); L(H)), i.e., M is bounded 
and analytic. 

Remark 2.3. The notion material law is motivated by several examples of mathematical 
physics, since it turns out that all material parameters, such as density, conductivity, per- 
meability etc. can be incorporated into the operator M(8q ) (see |16| for several examples). 
Thus, it is natural that M(d$ ' ) is forward causal. Since the operator M(8q ' ) is linear and 
it commutes with the translation operators Th, mapping u G H gt a(M.; H) to t \— > u(t + h) 
for h G K, causality can be characterized via the requirement that sptM(9^"^)n C]0, oo[ if 
spt u C]0, oo[, where by sptg we denote the support of a function g G L2 ) i oc (^; -ff)- This, how- 
ever, can be characterized by the analyticity and boundedness of the mapping M , employing 
a Paley- Wiener-type result (see e.g. |2Q^ Theorem 19.2]). Moreover, note that due to the 

boundedness of M, the operator M (Oq^ J becomes a bounded operator on H e! o(K;H). 



We are now able to state the solution theory for evolutionary equations. For the proof we 
refer to [T6j [23] . 

Theorem 2.4 (Solution Theory). Let A : D (A) C H — >• H be a maximal monotone linear 
operator. Moreover, let r > and M G T-L°°(Bc(r,r); L(H)) and assume that the solvability 
condition is satisfied: 

3c> Vz G B c (r, r) : <Kc z~ x M{z) > c. (4) 

Then for each g > i i/te problem (d|) is well-posed in H 8t o(B.; H), i.e. for each g > ^ f/ie 
operator 8q^M (8q X ) + A is boundedly invertible and has a dense range. Moreover, the solution 

operator ^8q^M{8q X ) + is forward causal. 

-i 



Remark 2.5. Note that the solution operator ^(?o it ,M(<9 *) + commutes with the time- 

derivative <9o,g- This yields, that for right hand sides / G D(8q ) for G N in |2| the cor- 
responding solution « also lies in D{8q ), i.e. the solution operator preserves "temporal" 
regularitjU. 



2 We denote by Bc(x, s) the open ball in C with center x £ C and radius s > 

3 Indeed, one can show that the solution operator (^do te M{d^ 1 B ) + extends to the whole Sobolev-chain 

(f/fc(<9o,g)) fegz associated with the derivative do, e , see |17j . which yields a solution theory for distributional 
right-hand sides. 
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3 An abstract condition for exponential stability 



In this section we show that under certain constraints on the function M, the corresponding 
evolutionary problem is exponentially stable. Although in the literature, exponential stability 
usually means that the solution of an initial value problem decays with an exponential rate 
as time tends to infinity, we like to introduce a slightly weaker notion within our framework, 
which however yields the desired decay, if the source term / is regular enough (see Remark 
13.21 (a)). Moreover, since the framework of evolutionary equations covers different types of 
equations, where initial values won't make sense, such as differential-algebraic or even pure 
algebraic equations or equations with memory effect, where a given pre-history would be 
more appropriate then an initial value, we cannot treat initial value problems within this 
general approach. However, in concrete examples we can reformulate initial value problems 
as evolutionary equations with a modified right hand side (see Subsection 4.1) such that the 
following results still apply. 

Definition 3.1. Let A : D(A) C H — > H be a maximal monotone linear operator and 
M G / H' X) {B£{r^r)]L{H)) for some r > 0, which satisfies the solvability condition (j4|). Let 

g > We call the solution operator ^3o i£ iM(9q"^) + A^j of (|2|) exponentially stable with 

stability rate v$ > 0, if for each < v < uq and / G H^ U) q(R; H) n H Qt o(R; H) the solution u 
of ([2|) satisfies 

u= (do^Mid^ + Ay 1 f £ P| H^o(R;H), 

— U<fJ,<Q 

which especially implies J R e 2>lt \u(t)\ dt < oo for all < fi < v. 
Remark 3.2. 

-i 



a) Let q > and ydo tS M(d ^) + Aj be exponentially stable with stability rate v$ > 0. 

Moreover, assume that / G H- u fi(M.; H) n H e fl(R;H) and / G D{8q )Q ) such that do jQ f £ 

-l 



H_ U)0 (R;H) n H efi (R;H) for some < v < u . Then u = ( d ^M{d^ Q ) + A) f also 

lies in D(<9o, e ) and do tQ u = (do^M^d^) + A^j <9o,£>/ (compare Remark I2.5p . By the 
assumed exponential stability, we get that e^ m u G L2(R;H) and e^ m do t gU G L2(M.;H) 
for each < fi < v. The latter yields e^ m u G W%(R;H). Indeed, for <p G C™(R;H) we 
compute 

(^ m u + e» m d , g u\<P) L2{R;H) 

(n(t)I^V(t)) dt + j ((9o i£) n)(t)| e ( 2 ^)V(t)e" 2e ') dt 

K 

(n(t)I^V(t)) dt + {do, g u\e^ m <P) HeMR . iH) 



(u(t)|^V(i)) dt + (u 
+ (u 2g U 2e+ ^ m 



d , e ( e^ + ^ m c 



H e , (R;H) 
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3 An abstract condition for exponential stability 



{u(t)\u,e^(f){t)) dt- J (2g + fi){u(t)\e^ +ll)t 4>(t))e- 2st dt 

R 

(n(t)|e (2 ^ t) V'(i))e" 2ei dt + J 2g(u(t)\e^ + ^ t (/)(t))e- 2et 

(«(t)|e( 2 ^V(t))e- 2 »* 
= -(e^ n |0') L2(M;H) . 

Thus, we obtain e^*|it(t)| — >• as t tends to infinity for each < fx < v due to Sobolev's 
embedding theorem (see e.g. [U p. 408]), i.e. u decays exponentially with a decay rate less 
than v. 

(b) If / G ff Ml o(K; H) n i^, (R; i?) for some g G R with //>£>, then / G H ufi (R; H) for all 
^ G [fiNA 4 ]- Indeed, we estimate 

oo 

J \f(t)\ 2 e~ 2vt dt= J \f(t)\ 2 e~ 2vt dt + J |/(t)| 2 e- 2 ^ dt 

R -oo 

oo 

< J \f(t)\ 2 e~ 2 ^ dt + J \f(t)\ 2 e- 2 ^ dt 



-oo 
- I/I# M>0 (R;H) + l/lif fl , (K;/r)" 

We now give conditions for the function M and show that they yield the well-posedness and 
exponential stability for the corresponding evolutionary problem. 

Hypotheses. Let i/q > 0. We assume that 

(a) M : C \ B c [- ^, ^] -> L(#) is analytic^, 

(b) for each r > and < v < z^o the function^ 

B C (r, r) \ {u" 1 } 3 z i-> (1 - vz) M (z(l - vz)~ l ) 

has a bounded and analytic extension to Bc(r,r), 

(c) for every < v < uq there exists c > such that for all z G C\ Be \— 

d\tz~ 1 M{z) > c. 

Remark 3.3. Let M satisfy the assumptions above and let r > 0. Then the restriction of M 
to Bc(r,r) is an element ot 71°° (Bc(r,r); L(H)). Indeed, the analyticity is clear from (a) and 
the boundedness follows from (b). Hence, together with (c), it follows that the problem (|2|) is 
well-posed for every g > according to Theorem 12.41 

4 We denote by Bc[x, s] the closed ball in C with center i£C and radius s > 0. 
J Here and further on we set B c (r,r) \ {(T 1 } := B c (r,r). 



10 



We now state some auxiliary results, which in particular imply that the solution operator 
\ _1 

9q,qM ypog^j + A j f° r an evolutionary problem does not depend on the particular choice 



of q. These results can also be found in |17| p. 429 f.]. However, for sake of completeness we 
present them again with a slightly modified proof. 

Lemma 3.4. Let g,fj, G R with fj, > g and set U := {z G C | 9\t z G [g, [/]}. Moreover, let 
f : U H be continuous on U and analytic in the interior ofU, such that f^" f(iR+s) ds — > 
as R — > ±oo. Then 

R R 

limsup / /(it + fi) dt = limsup / f(it + g)dt. 
R— >oo J R—>oo J 

-R -R 

Proof. According to Cauchy's integral theorem we have 

R » I R A* \ 

i J /(it + g) dt + J f(iR + s) ds - I i J /(it + dt + J /(-LR + s) ds J =0 
—R e \-R g I 

for each R > 0. The assertion now follows by taking the limes superior as R tend to oo. □ 

The next lemma shows that we can approximate a function, which belongs to two different 
exponentially weighted L2-spaces by the same sequence of test functions with respect to both 
topologies. 

Lemma 3.5. Let ft/j£M and f G H e> o(M.; H) n -£f Mi o(R; H). Then, for each e > there exists 
4> G C™(R;H) such that 

max{| / - (/)\ He , (R;H), \f ~ 0\h^ o (R;H)} < £■ 
Proof. Let e > 0. Then we choose N G N such that fpf := fX[-N,N] satisfies 
max{|/ - fN\H e , (U;H), 1/ - /iv|i? M>0 (R;i?)} < ^ 

We denote by (V>fc)fceN G C C °°(R) N the Friedrichs-mollifier (see e.g. (TDJ Chapter C.4]). Then, 
for each k G N we have that spt ipk * /at Q [~ N — 1, N + 1]. Now, we choose feo large enough 
such that 

\„i. .f f I ^ -2max{|£|,|«|}(7V+i) ' 



Then the function ijj^ * /jv G C£°(K; i?) has the desired property. □ 

Lemma 3.6. Let £ R with \i > g and set U := {z £ C 9\e z G [g, //]}. Moreover, let 
f G H gfi (.^;H) n H^ {R;H) and T G H°°(U;L(H)) n C b (U;L(H)) (i.e. T is bounded and 
continuous on U and analytic in the interior ofU). Then 

(£*T(im + (t) = (£*T(im + (t) 

/or almost every t G R. 
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3 An abstract condition for exponential stability 

Proof. According to Lemma 13.51 it suffices to prove the assertion for test functions. So let 
cj) G C C °°(]R; H). We show that the function 

Ubz\-¥ e zt T{z) ((£?r c ^) (3m z)) G H 

satisfies the assumptions of Lemma f3.4l Indeed it is continuous in U and analytic in the interior 
of U as a composition of analytic functions. Furthermore, we estimate for s G [g,(i],(£l: 

\£ S J>(0\ < A= J e~ rs \Hr)\ d ^ "7=/ 1^)1 dr, 



where C := sup{e~ rs \s G [g, fj],r G spt0}, which shows that the function 

U 3 (C<R tz 4>) (3m z) 

is bounded. Moreover, due to the Riemann-Lebesgue lemma, we get that (C s <p) (R) — > as 
R =too for every s G R. Therefore, according to Lebesgue's dominated convergence theorem, 
we deduce that 

J e S lR+s ^T(iR + s) (C s (f)) (R) ds < max{e et , e^T^ j \(C s <p)(R)\ ds -> (12 ±oo). 

Thus, by Lemma [3.4| we get that 

J e^ s+ ^T(is + g) (L e <t>) (s) ds = J e^ + ^T(is + n)L^(s) ds, 

R R 

which yields the assertion. □ 
We are now able to prove our main theorem. 

Theorem 3.7. Let A : D(A) C H —> H be a maximal monotone linear operator and M a 

mapping satisfying the hypotheses above for some z^o > 0. Then for each g > the solution 

-l 



operator ydo t gM(d ^) + Aj is exponentially stable with stability rate Uq. 
Proof. Let g > 0, < v < v and take / G H- vfi (R; H) n H efl (M.- H). We set 



u := 



do, e M(d^) + A) f 

and we have to show that u G ff«o(R; H) for each Let < rj < g + v. We define 

iV» = (1 - vz)M (z (1 - vzY 1 ^ 

for z G i?c (^>2n) \ i^ 1 }- Note that according to hypotheses (b), TV has an extension 
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N G H°° [Be l± ± ; L(H) . Moreover 



<Re z~ x N(z) = <Rt z~ l N{z) = IHe (1 M \T~^) - c ^ 



for every z G Be ( ^ j \ {v an d some suitable c > 0, since z (1 — vz) 1 ^ 2?c [— ^tj, ^tj] • 
Due to the continuity of N, estimate ([5]) also holds for z = Thus, according to Theorem 



2.41 we obtain a solution 



« := (5 0ir) iV(5 -J) + e um f G iT^R; ff), 

where we have used that e^ m / G ff ,o(R;-£0 H ff e+I/)0 (M; -ff) C H vfi (R;H) (see Remark 
(b)). We apply Lemma [3761 to T(z) = (zN (z~ l ) + A)~ l for z G C with <Kez > r/ and get that 



„ = £; (( im + ,)Ar(_L_) + ^ 



-1 

r -'l 



C v e vm f 



£*+» ( im + V + ")N[ . - - + A Cn+ueTtf. 

im + Tf + 1/' 



-1 



Using (it + 77 + v) N (n4+7j) = (it + r/)M for t G K, we derive 



Again, applying Lemma 13.61 to T(z) = (zM(z 1 ) + A) for z G C with 9^ez > minj^, 77} 
we get that 



<T vm v = Cl \ (im + r?) M ( — L- J + A J £,/ 



im + 77 



-1 



(-^) 



which gives u G H„_j,(R; -ff). Since 77 G]0, g + z^ was chosen arbitrarily, we derive the assertion. 

□ 



4 Examples 

In this section we illustrate our results of the previous section by three different types of 
differential equations, which, however, are all covered by the abstract notion of evolutionary 
equations. We emphasize that we do not claim that in the forthcoming examples the stability 
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4 Examples 



rates are optimal under the given constraints nor that an exponential decay could not be 
obtained under lesser constraints. But we emphasize that our approach provides a unified 
way to study exponential stability of a broad class of differential equations. 
We begin to study a class of differential-algebraic equations, where the material law is of the 
simplest form. Moreover, we provide a strategy how to deal with initial values problems for 
this class. As a second example we treat a partial differential- algebraic equation with finite 
delay. We conclude this section with an example of a parabolic integro-differential equation 
with an operator-valued kernel. 



4.1 Differential-algebraic equations of mixed type 

It turns out that in applications, the material law is often of the form M(9q~*) = Mq + Oq^Mi 
(see for instance |16 | I17 | [2T]). where Mq, M\ € L(H). The corresponding evolutionary equation 
is then of the form 

(d , e M + M l + A)u = /, (6) 

where A : D(A) C H — > H is a maximal monotone linear operator. In order to obtain the 
well-posedness of this evolutionary equation we require that Mo is selfadjoint and strictly 
positive definite on its range, while £Re Mi is strictly positive definite on the null-space of Mo 
(see |16| EH E3] for the proof of well-posedness). In order to obtain an exponential stability 
for this problem, we require that M\ is strictly positive definite on the whole space H. 

Theorem 4.1. Let A : D{A) C H — >• H be a maximal monotone linear operator and Mo, Mi £ 
L(H) such that Mq is selfadjoint and strictly positive definite on its range, and £Re Mi > c > 0. 
Then for each g > the solution operator (<9o l£ Mo + Mi + A) is exponentially stable with 
stability rate pfgji- 

Proof. We have to verify the hypotheses (a)-(c) for the function 

M{z)=M + zM 1 (z€C). 

Obviously the assumption (a) holds. Let now r > and v > 0. Then we compute 

(1 - vz)M (z(l - vz)~ l ) = (1 - vz)M + zM x 

for z G Bc(r, r) \ {u^ 1 }, which shows (b). Let now v g]0, pf^jft- I n order to show (c) we note, 
that for z£C \ P>£ [ — there exists t € M and g > —v such that z _1 = it + g. Thus, for 
g > we can estimate 

?te z~ 1 M(z) = qM + 9\e Mi > c, 
while for g G ] — v, 0[ we estimate 

D\e z~ 1 M(z) = gM + <Ke Mi > -z/||M || + c> 0. 

Thus, the assertion follows by Theorem 13.71 □ 

To illustrate the versatility of our approach, we discuss the following simple example of an 
evolutionary equation, whose type (elliptic, parabolic or hyperbolic) changes on different parts 
of the underlying domain. A similar example was also discussed in |18| . 
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4.1 Differential-algebraic equations of mixed type 



Example 4.2. Let f2 C M n and f&o>^i C SI be measurable disjoint subsets with positive 
Lebesgue-measure. We consider the evolutionary equation 

*-( XB, : xB, «) + (i!2) + U*-))(:)-(;)- (7) 

where c > 0. The differential operator div c is defined as the closure of the operator 

div|cco (tl )» : c?(n) n C L 2 (n) n L 2 (fi) 



)ie{l,...,n} ' 



i=l 



where we denote by c\ the partial derivative with respect to the z-th coordinate. The operator 
grad is defined as the negative adjoint of div c , i.e. 

grad := - (div c )* . 

This operator is just the usual weak gradient on L^iP) with domain H l {Q). Note that then the 
operator matrix ( ^ J is a skew-selfadjoint operator (and hence maximal monotone) 



grad 

on H := L/2(fl) © L2{£l) n . Moreover, the operators 

V Xn J \0 c 

satisfy the assumptions of Theorem 14. 1 1 and hence, the solution operator is exponentially stable 
with stability rate c. 

Although this example seems to be quite easy, it seems hard to attack the problem of solving 
([6]) by using semigroup techniques. The reason for that is that ([6]) changes its type on different 
parts of the domain f2. Indeed, on we obtain a hyperbolic problem of the form 



1/ \0cJ \ grad J J \q J \g 
while on S}\ the problem becomes parabolic, namely 

(l o) + (o c) + (grad d C )) (g) = {{ 
which yields, in case of g = the parabolic differential equation 

do, g v - c 2 div c grad v = f. 
On the remaining part Q \ (CIq U Oi) the problem is elliptic 

c0\ / div c \\ /t;\ = // 

o c y + Vs rad yy \q) \a 

Note that we can treat this problem, without requiring any explicit transmission conditions on 
the interfaces dilo and d£l\ and without imposing regularity assumptions on these boundaries. 
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Now, we present a possible way to tackle initial value problems for equations of the form ([6]). 
Consider the following initial value problerrjf] 

(d Q , e M + Mi + A) u = f on]0,oo[ (8) 
M u(0+) = MqUq 

for M , Mi, A as before, / G iT_„ i0 (K; ff) n ff e ,o(K; #) for some v, g > with spt / C [0, oo[ 
and uq G -D(^4)- One way to deal with such problems is to consider the evolutionary equation 



(do )£? Mo + Mi + A) v = f - X[o,oo[( m ) MlU ° ~ X{o,oo{( m ") Au o 



on 



for the new unknown v := u — X[o,oo[{ m ) u o given by v(t) = u(t) — X[o,oc[(t) u o f° r almost every 
t £ R. Then the right-hand side belongs to H Sj q(R;H) for positive g, but does not decay 
if / decays. Hence, this approach can be used for the issue of well-posedness but it is not 
appropriate for exponential stability. 

Instead, we consider an alternative problem for the unknown 

v := u — 4>(m)uo, 

where <j) is given by 



l ifte[0,l], 
2-t if*€]l,2[, 
otherwise. 



It is clear that, if u satisfies (JH]), then v satisfies the equation 

(d 0>g M + Mi + A) v = f + X]i,2[(m)M u - ^(m)M lUo - (j)(m)Au -■ g (9) 

and vice versa. Since the function X}i,2[( m )MoUo— (f)(rn)M\UQ— 4>(m)Auo belongs to iJ^ofM; H) 
for every \i G R, we obtain g G H gfi (R;H) n H- Vj0 (R;H). Thus, Theorem O applies to © 
and we get that (|9|) is well-posed and v G Hno(M;H) for each fj, g] — v, g]. This gives that 
u G Hnfi(M.; H) for each |iE]- v, q], since <fr(m)uo G (~) ^H^oiM.; H). It is left to show in 
which sense Mqu attains the initial value MqUq. By ([9]) we get that 

d 0>e M v = g - Miv - Av 

and the right hand side belongs to H e o(M.; H-i(A + l)), where H_\(A+1) is the extrapolation 
spaced associated with A + 1. Thus, a version of Sobolev's embedding theorem ( )17[ Lemma 
3.1.59], |131 Lemma 5.2]) yields that Mqv is continuous as a function which attains values in 
iZ_l(A+l). Furthermore, due to the causality of the solution operator (9o i£ Mo + Mi + ^4) , 
v is supported on [0, oo[, since spt g C [0, oof. This yields Mqv(0+) = Mqv(0—) = and hence, 
since 0(0+) = 1 we get that Mqw(0+) = MqUq, where the equality holds in H_\(A + 1). 



6 Note that it only makes sense to prescribe an initial value for MoU yielding an initial value for the part of 
u lying in the range of Mo. 

7 For a boundedly invertible linear operator C : D(C) C H H the extrapolation space H-i(C) is given as 
the completion of H with respect to the norm | • \h_ 1 (c) defined as |a;|jj_ 1 (c) : = IC" 1 ^;!^ for x € H (see 
[171 Section 2.1]). 
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4.2 Linear partial differential equations with finite delay 



4.2 Linear partial differential equations with finite delay 

As a second example we study a differential equation with finite delay of the form 

(do, e M + T h + M l + A) u = /, (10) 

where Mo, Mi G L(H) such that Mq is selfadjoint and non-negative, ^fKzM\ > c > 1, A : 
D(A) C H — > H is linear and maximal monotone and is the translation operator with 
respect to time, i.e. (r^it) (£) = uit + h) for £ G R and some /i < 0. We will prove that under 
these assumptions, the corresponding solution operator is exponentially stable and we give an 
estimate for the stability rate. A similar problem is treated in [2J Example 4.14], where the 
well-posedness is shown via semigroups and a criterion for the exponential stability is given, 
using the spectral mapping theorem for eventually norm continuous semigroups (cf. [9J p. 
280]). 

Before we state our stability result for (fTU|) . we need to inspect the operator a bit closer. 
Lemma 4.3. Let We define the operator 

T h : H ei0 (R;H)^H et0 (R;H) 

u i — y (£ i — y u(t + h)) . 

Then r h G L(H g>0 (R; H)) with \\T h \\ = e eh and r h = C* Q e^ m+ ^ h C e . 



Proof. Obviously, Th defines a bounded linear operator on H s< q(R; H). For (j) G C£°(R;H) we 
compute 



C e (r^) (t) = -j=f e- ist e-e s 0(s + h) ds 



e {it+ e )h_^ = l e -bt e -o»j > ( 8 \ d s 



'2ix 
R 

= e^ h (C g c/>) (£) 

for each £ G M, which gives = C*e^ m+e ^ h C g . Moreover, by this unitary equivalence we get 
that 

\\r h \\ = ||e (im+ ^|| = e eh . □ 



Using this Lemma, we are able to write (|10p as an evolutionary equation with 

M(z) = M + ze z ~ lh + zM 1 , (11) 

which is clearly analytic and bounded on balls of the form B<c(r,r) for r > 0, if we require 
that h < 0. This restriction is natural, since for h < the operator is forward causal, while 
for h > it is backward causal. 

Theorem 4.4. Let A : D(A) C H — >■ H be a maximal monotone linear operator, Mo, Mi G 
L(H) such that Mq is selfadjoint and non-negative and D\eMi > c > 1. Moreover let h < 0. 
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Then for each g > the solution operator [Bq^Mq + r/j + M\ + ^4) 1 is exponentially stable 
with stability rate vq > satisfying 

u \\M \\ + e~ Uoh = c. 



Proof. We have to show that M given by (Hip satisfies the hypotheses (a)-(c) of Section 3. 
Obviously M is analytic on C \ {0}, which shows (a). Let now r > and v > 0. As 

(1 - vz)M [z{\ - vz)~ x ) = (1 - uz)M + ze^' 1 '^ + zM x 

for z G Bc(r,r) \ {z^" 1 }, we see that z i— >■ (1 — vz)M {z(\ — vz)~ 1 } has an analytic extension 
to .Be (r, r) . Moreover we estimate 

sup le^ -1 -"^! = sup e^ h , 

zeB C (r,r) g> ±. 

which is finite, since h < 0. Thus, hypothesis (b) is also satisfied. For showing (c), let v G]0, v§\ 
and z G C\Bc [—^7, • Then there exists £ > —1/ and t£R such that z" 1 = it + £>. If £> > 
we estimate 

<He 2 _1 M(;y) = qM + e e/l cos(t/i) + V\tM x > c - 1 > 0, 
and for the case g g] — v, 0[ we get 

We3f -1 M(z) = pM + e eh cos(t/i) + <HeMi > -i/||M || - e~ vh + c> 0, 
since v < z^. This proves (c) and thus, the assertion follows by Theorem 13.71 □ 

Remark 4.5. Note that M$ in (jlOp is allowed to have a non-trivial null-space, which could 
also depend on the spatial variable (compare Subsection 4.1). Thus, Theorem 14.41 also covers 
a certain class of differential-algebraic equations with delay. 



4.3 Parabolic integro-differential equations 

We consider the following parabolic integro-differential equation 

d Q , e u + Bu- C * Bu = f, (12) 

where B : D{B) C H — >■ H is linear such that ^4 := B — c is maximal monotone for some 
c > 0, C : [0,oo[— )• L(H) is a weakly measurable function such that i h-> ||C(t)|| is measurable 
and there exists i/ > with / °° ||C(t)||e i/ot dt < 1. We set U := {z G C 3m z < z/ } and 
define the complex Fourier-transform of C by 

00 

e- itz C(t) dt (z G U), 
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4.3 Parabolic integro-differential equations 



where the integral is meant in the weak sense. Note that C : U — > L{H) is continuous and 
bounded on U and analytic in the interior of C7|l. The well-posedness and asymptotic behavior 
for equations of the form (|12|) . including non-linear perturbations, were discussed in several 
works (e.g. [5j [Ml [22] and [6l [191 El El [22] for a hyperbolic version of the problem), imposing 
additional constraints on the kernel C. 

Following [22 j we are led to assume that C satisfies the following conditions: 

1. C(t) is self adjoint for almost every t G K, 

2. C{t) and C(s) commute for almost every t, s G R, 

3. for all t G M we have 

t3md(t + iv ) <0. (13) 



Before we can state a stability result for problems of the form (|12p . we recall some properties 
of the convolution operator C * . 

Lemma 4.6. We denote by S(M; H) the space of simple H -valued functions. Then for q > — uq 
the operator 

C* : 5(R; H) C H efi (R; H) -> H e>0 (R; H) 



u i — y j 1 1 — ^ / C(t - s)«(s) ds 
\ R / 



is bounded and linear with \\C * \\l(H (R-,H)) — Jo° \\C(t)\\e' yot dt and can therefore be extended 
to H e o(R; H). Moreover, for u G H e o(R;H), q > — vq we have 



{C e {C * u)) (t) = V2vrC(t - ig) (£ s ii(t)) (14) 

for almost every t £ 1. 

Proof. A proof for the first assertion can be found in |221 Lemma 3.1]. The proof of formula 

(|14p is straight forward and we omit it. □ 



According to Lemma [4.61 the operator (1 — C*) is boundedly invertible in iJ gi o(M; H) for each 
Q ^ — vq. Therefore, instead of considering ()12p we can study 



(d Q>e (l - C*)- 1 +B)u = (l- C*)" 1 /, (15) 

or equivalently 

(do, e (l - C*)- 1 + c + A) u = (1 - C*)- 1 /. 
Note, that this as an evolutionary equation of the form ([2]) where M is defined by 



M(z) = (1 - v^CKz- 1 ))- 1 + cz, (z€C\ B C (-^-, 



(16) 



3 Here we use the fact that scalar analyticity and local boundedness on a norming set yields analyticity (see 
[T2l Theorem 3.10.1]). 
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where we have used Lemma 14.61 The next Lemma shows that (|13ft already implies that the 
same condition holds if one replaces — vq by g for arbitrary g > —uq. 

Lemma 4.7. Assume that C satisfies the conditions 1., 2. and 3. above. Then for every 
Q > — we have 

1 3m C(t - ig) < 0. 



Proof. The proof can be done analogously to the one of |22| Lemma 3.7]. 



□ 



We now state our stability result for (1151) . 

Theorem 4.8. Let A : D (A) C H — > H be maximal monotone and linear and let c > 0. 
Moreover, let C : [0, oof— > L(H) be weakly measurable, such that t i— > ||C(£)|| is measurable 
and there exists uq > such that \\C(t)\\e uot < 1 and C satisfies the conditions 1., 2. and 



3. from above. Then for each g > the solution operator (^<9o,e(l — C*) 1 + c + A 
and is exponentially stable with a stability rate u\ G]0, uq] satisfying 



exists 



u x 1 



\\C{s)\\e^ ds 



< c. 



Proof. Let v\ g]0, v Q ] such that v x (l - J °° ||C(s)||e 1/lS ds) 1 < c. We prove that M given by 
(|16p satisfies the hypotheses of Section 3. The assumption (a) is clear, since ||v / 2vrC'(— iz" 1 )!! < 
1 for each z G C\Bq — ^ . Let now r > and < v < V\. Then for z G i?c(r, r)\{^~ 1 } 
we compute 



(1 - vz)M(z(l - vz)~ l ) = (1 - vz) ( (l - V2^C {-iz~ l {l - uz))) + cz{\ - vz)~ l 



(1-vz) [1 - V2^C (-{(z- 1 - v))) + 



-i 



cz, 



which has a holomorphic extension in v . Noting that for each z G Bc{r,r) we have that 



it + £ for some g > ^ , t G M we estimate 



oo 



-(e-f)s 



||C(s)|| ds 



< / e I/0S ||C(s)|| ds < 1. 
o 

Hence, the extension of (1 — vz)M(z(l — vz) -1 ) to Bc(r,r) is indeed bounded for each r > 
0, v G [0, i^i [. We now show that M satisfies the assumption (c) on C \ Bq 2lT ■ We 

follow the strategy of the proof of |22| Lemma 3.8]. Let z G C \ -Be — 2?a ' 2^1 ' Note that 
then there exists g > — V\ and t G K such that z" 1 = if + g. We set L> := |1 — y/27rC(t — i^)!" 1 , 
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4.3 Parabolic integro-differential equations 



which is well-defined according to Lemma [4.71 Note that 

(1 - V2^C(t - ig))" 1 = 25 2 (1 - V2^C(-t - ig)), 
where we have used assumption 1. Moreover, due to assumption 2., we get that 

D 2 (l - V2nC(-t - ig)) = 25(1 - V2nC(-t - ig))D. 
Thus, we obtain for x £ H 

9\e(z- 1 M(z)x\x) = JRe / z' 1 ((1 - {-iz' 1 ))' 1 + cz) x 



x 



<Ke ( ( g d\t ( 1 - V2^C(-t - ig)) + y/2nt 3m C(-t - ig) ) Dx 



Dx ) + c\x\ 



> g <Re / D ( 1 - v / 2^C(-t - ig)) Dx 



xj + c\x\ . 

If g is non-negative, the latter term can be estimated by c. For negative g we observe that 

-i 



25<Ke (1 - V2nC(-t-iQ)) 25 





< 




< 




~ ] 




<-\ 







1 - V2^C(-t - ig) 
1 



1 - ||V27rC(-t-if 



-i 



and hence, 



g <Re ^ 25 (l - V2^C(-t - igf) Dx 



x) + c\x\' 2 > [ g ( 1 - J e" ls ||C(s)|| ds | + c | |x| 

> | -z/x I l-Je^ s \\C(s)\\ dsj +cj |x| 2 , 

which shows that M satisfies hypothesis (c), according to the choice of v\. Thus, the assertion 
follows by Theorem 13.71 □ 



Remark 4.9. Theorem 14.81 gives the exponential stability for equation (|15p . This also yields 
the exponential stability of the original problem (|12p . since the operator (1 — C*)" 1 leaves the 
space 22-^0 (R; -22) for all v < uq invariant. Indeed, observing that 

e um (l - C*)" 1 = (1 - {e um C) *) _1 e um 

we obtain for / G 22_ I/>0 (R; H) n 22" ei0 (R; 22") 

e ™(l _ C7*)" 1 / = (1 - ( e ^ m C) *) _1 e um f G L 2 (M; 22). 
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